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9 Abstract. The nonlinear turbulent interactions between acoustic grav- 

10 ity waves are investigated using two dimensional nonlinear fluid simulations, 
u The acoustic gravity waves consist of velocity and density perturbations and 
12 propagate across the density gradients in the vertical direction in the Earth's 
o atmosphere. We find that the coupled two component model exhibits gen- 

« eration of large scale velocity potential flows along the vertical direction, while 

is the density perturbations relax towards an isotropic random distribution. The 

16 characteristic turbulent spectrum associated with the system has a Kolmogorov- 

17 like feature and tends to relax towards a k~ 5 ^ 3 spectrum, where k is a typ- 

18 ical wavenumber. The cross field diffusion associated with the velocity po- 

19 tential grows linearly and saturates in the nonlinear phase. 
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1. Introduction 

20 Studies of acoustic gravity waves (AGWs) [Yeh & Liu , 1974; Frits & Alexander , 2003; 

21 Mckenzie & Axford , 2000; Stenflo , 1987; Jovanovic et al. , 2002] in the Earth's middle 

22 atmosphere, in the solar atmosphere, as well as in planetary magnetospheres have been 

23 motivated by the need to obtain accurate predictions of the dynamics of the atmosphere 

24 under various meteorological conditions, and include the profiles of the density, pressure, 

2 5 and the presence of shear flows in the winds. The AGWs are low-frequency disturbances 

26 associated with the density and velocity perturbations of the atmospheric fluid in the 

27 presence of the equilibrium pressure gradient that is maintained by the gravity force. The 

28 frequency ui of the AGWs is given by the dispersion relation [Stenflo , 1987; Jovanovic et 

2 9 al. , 2002; Stenflo , 1996] 



w 2 = k 2 y g /(ki + e z + i/4H*), (i) 

where k x and k z are the components of the wave vector along the x and the z axis in 
a Cartesian coordinated system, H is the density scale height, and the squared Brunt- 
Vaisala frequency is 

.j. _ ( dPo 1 , 1 d p \ ( d Po\ ( 2 ) 
9 \ dz 7p Po dz j \dz J ' 

30 where po(z) and po(z) are the equilibrium density and pressure that are inhomogeneous 

31 in the vertical direction (denoted by the coordinate z), and 7 « 1.4 is the adiabatic index. 

32 At equilibrium, we have dp /dz = —pog, where g = — gz is the acceleration force on the 

33 atmospheric fluid due to gravity, and z is the unit vector in the vertical direction. If 

34 oo 2 g < 0, it turns out from (1) that the AGWs are unstable. When the amplitudes of the 
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AGWs become large, the modes start interacting among themselves. The mode couplings 
are governed by [Stenflo , 1996] 



D 



dx 



(3) 



Dt 



dx 



and 




(4) 



s? where ip(x,z) is the stream function, x( x , z ) is the normalized density perturbation, 

as D/ Dt = (d/dt) + v • V, v = z(dip/dx) — x(dip/dz) is the fluid velocity, and V 2 = 

so (d 2 /dx 2 ) + (d 2 /dz 2 ). Time and space coordinates in Eqs. (3) & (4) are normalized re- 

* spectively by their typical values t and £ . The other variables are normalized as follows; 



« ?/>M, = f,Wo/k = X',t 2 uj g = oj' g ,H/£ = H',£ V = V,t D/Dt = D/Dt'. The 



42 primes have been removed from Eqs. (3) & (4). Stationary nonlinear solutions of (3) 

43 and (4) in the form of a double vortex and a vortex chain have been presented by many 

44 authors [Stenflo , 1987; Jovanovic et al. , 2002; Stenflo , 1996; Stepanyants , 1989; Heijst 

45 & Kloosterziel , 1989]. 

The coupled nonlinear equations conserve total energy E in the absence of sources and 
sinks. In the normalized units, it is 



46 It is noteworthy from the linear dispersion relation that the group and phase speeds of 

47 the AGWs are negligibly small when k x 3> k z . In such cases the fluctuations merely 

48 oscillate with the Brunt- Vaisala frequency. In the other limit k x k z , the AGWs are 

49 dispersive and have anisotropic propagation. Equations (3) and (4) possess the nonlinear 
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so terms e y x VV 2 ^ • Vx and e y x • V% that are analogous to the polarization and 

si diamagnetic nonlinearities in inhomogeneous drift wave turbulence. The former cascades 

52 energy towards larger scales, while the latter leads to the formation of smaller scales due 

53 to forward cascades. In the next section, we shall explore the nonlinear interactions of 

54 acoustic gravity waves. 

55 The nonlinear interaction between finite amplitude AGWs in the Earth's atmosphere 

56 [Gill , 1982] leads to energy transfer, resonantly and non-resonantly between the high- 

57 and low-frequency parts of the spectrum. Resonant interactions between the AGWs play 

58 also a key role in a rotating atmosphere [Axelsson et al , 1996]. Many analytical theo- 

59 ries [Stenflo , 1994] describing nonlinear couplings between different wavelength acoustic 

60 gravity modes predict the formation of vortices, which can be important for weather pre- 

61 dictions. Gardner Gardner [1994] explains vertical wave number spectrum by means of a 

62 scale-independent diffusivity by assuming the damping effects of molecular viscosity, tur- 

63 bulence, and off-resonance wave-wave interactions. Similarly, [Beres et al , 2004] presented 

64 an implementation of a physically based gravity wave source spectrum parameterization 

65 over convection in Global Circulation Model. 

66 With the objective of developing a self-consistent turbulent spectrum of AGWs we, in 
this paper, present the turbulent properties of nonlinearly interacting AGWs that are 

68 governed by Eqs. (3) and (4). Specifically, we investigate by means of computer simula- 

69 tions the statistical properties of the turbulent spectrum arising from dual cascading. The 

70 fluid diffusion across the density and pressure inhomogeneity directions in the presence of 
?! coherent structures is also deduced. Furthermore, we consider the statistical turbulence 
72 properties of the nonlinearly interacting AGWs. Specifically, we perform computer sim- 
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73 ulation studies of the acoustic gravity mode cascading as well as the resulting structures 

74 and turbulent spectra. The diffusion of the fluid across the density and pressure gradients 

75 in the presence of large scale acoustic gravity mode structures is also investigated. 

2. Nonlinear simulations 

76 To investigate the nonlinear mode coupling interaction and turbulence aspects of acous- 

77 tic gravity waves, we have developed a two dimensional spectral code to numerically in- 

78 tegrate equations (3) and (4). The 2D simulations are not only computationally simpler 

79 and less expensive compared with the full 3D simulations, but they offer significantly 
so higher resolution even on moderately sized small-cluster machines like the IGPP (Insti- 

81 tute of Geophysics and Planetary Physics, at University of California Riverside) Beowulf. 

82 The spatial discretization in our code uses a discrete Fourier representation of the turbu- 

83 lent fluctuations. The evolution variables use periodic boundary conditions. The initial 

84 isotropic turbulent spectrum is chosen close to a k~ 2 spectrum with random phases in all 
as the directions. The choice of such (or even a flatter than -2) spectrum treats the turbulent 
ae fluctuations on an equal footing and avoids any influence on the dynamical evolution that 

87 may be due to initial spectral non-symmetry. The equations are advanced in time using 

88 a Runge-Kutta 4 th order scheme. The code is made stable by a proper de-aliasing of 
as spurious Fourier modes and by choosing a relatively small time step in the simulations. 

90 Our code is massively parallelized using Message Passing Interface (MPI) libraries to facil- 

91 itate higher resolution in a 2D computational box. The computational domain comprises 

92 512 2 modes in two dimensions. Other simulation parameters in the normalized units, as 

93 described above, are H = 10, uj g = 0.02, and L = 4 x 4 (two dimensional box size). 
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94 All the fluctuations in our simulations are initialized isotropically with random phases 

95 and amplitudes in Fourier space. The initial fluctuations do not possess any flows or mean 

96 fields. The latter can be generated in our simulations by means of nonlinear interactions. 

97 Fourier spectral methods are numerically almost non dissipative compared to the existing 

98 finite difference methods and are thus remarkably successful in describing turbulent flows 

99 in a variety of plasma and hydrodynamic fluids. They also provide an accurate represen- 

100 tation of the fluid fluctuations in Fourier space. Because of the latter, nonlinear mode 

101 coupling interactions preserve ideal rugged invariants of fluid flows, unlike finite difference 

102 or finite volume methods. The conservation of the ideal invariants such as energy, enstro- 

103 phy, magnetic potential, helicity etc. in the turbulence is an extremely important feature 

104 because these quantities describe the cascade of energy in the inertial regime, where tur- 
ns, bulence is, in principle, free from large-scale forcing as well as small scale dissipation. We 
we include small scale dissipation, however, to distinctly extend the inertial range that helps 

107 push the spectral cascades further down to the smallest scales. The latter leads to an 

108 improved statistical average of spectral indices. The numerical validity in the simulations 

109 is checked by monitoring the integrable form of the energy of the system. 

no The initially randomly propagating acoustic gravity waves begin to interact linearly 

in during the early phase of the simulations. As the waves acquire larger amplitudes, they 

U2 begin to interact nonlinearly. During the nonlinear interaction phase, various eddies 

to mutually interact and transfer energy between the modes. The evolution of the modes 

n4 for density fluctuations and velocity potential fields are shown in Figs. (1) and (2) in 

us the (x, z) plane. It appears that the density perturbations (i.e. x ) have a tendency 

no to generate smaller length-scale structures, while the velocity potential cascades towards 
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U7 larger scales. This is consistent with the corresponding nonlinear terms in drift wave 

us turbulence. The co-existence of the small and larger scale structures in the turbulence is 

U9 an ubiquitous feature of various 2D turbulence systems. For example, in 2D drift wave or 

120 CHM turbulence, the plasma fluid admits two inviscid invariants, namely the energy and 

m the mean squared vorticity (i.e. irrotational velocity fields). The two invariants, under 

122 the action of an external forcing, cascade simultaneously in turbulence, thereby leading 

123 to dual cascade phenomena. 

124 In these processes, the energy cascades towards longer length-scales, while the fluid 

125 vorticity transfers spectral power towards shorter length-scales. Usually a dual cascade is 

126 observed in a driven turbulence simulation, in which certain modes are excited externally 

127 through random turbulent forces in spectral space. The randomly excited Fourier modes 

128 transfer the spectral energy by conserving the constants of motion in k space. On the 

129 other hand, in freely decaying turbulence, the energy contained in the large-scale eddies is 
wo transferred to smaller scales leading to a statistically stationary inertial regime associated 

131 with the forward cascades of one of the invariants. Decaying turbulence often leads to 

132 the formation of coherent structures as the turbulence relaxes, thus making the nonlin- 

133 ear interactions rather inefficient when they are saturated. It seems that the acoustic 

134 gravity waves exhibit turbulence in a manner similar to two dimensional hydrodynamic 
us turbulence in which the density and velocity potential fluctuations are analogous to the 
no hydrodynamic fluid energy and vorticity. The evolution of the turbulent energy associated 
137 with the modes (density and potential field) are shown in Fig. (3) where the linear and 
us nonlinear phases of the evolution are clearly marked by the curve. After the nonlinear 
139 interactions are saturated, the energy in the turbulence does not grow and remains nearly 
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«o unchanged throughout the simulations. Correspondingly, the energy transfer rate shows a 
i4i significant growth during the linear and initial nonlinear phases. However when the non- 
142 linear interactions saturate, the nonlinear transfer of energy in the spectral space amongst 
us various turbulent modes becomes inefficient and the energy transfer per unit time tends 

144 to become negligibly small as shown in Fig. (4). The power spectrum associated with 

145 the acoustic gravity wave turbulence exhibits a spectral slope close to —5/3 as shown in 
He Fig. (5). This is indicative of the turbulent eddies being the dominant processes in the 
147 spectral transfer and is consistent with the Kolmogorov-like phenomenology. 

Ms We finally estimate the turbulent transport coefficient associated with a self-consistent 

149 evolution of the small and large scale turbulent fluctuations. An effective electron dif- 

150 fusion coefficient caused by the momentum transfer can be calculated from D e ff = 

151 / °°(v(r, t) ■ v(r, t + t'))dt', where v is the fluid velocity, the angular bracket represents 

152 spatial averages and the ensemble averages are normalized to unit mass. Since the 2D 

153 structures are confined to a x — z plane, the effective diffusion coefficient, D e ff, essen- 

154 tially relates the diffusion processes associated with the translational motion of the fluid 

155 elements of our nonlinearly coupled system. We compute -D e // in our simulations from 
we both the velocity field fluctuations to measure the turbulent transport that is associated 
is? with the large scale flows. It is observed that the effective cross-field transport is lower 
us initially, when the field perturbations are Gaussian. On the other hand, the diffusion 

159 increases rapidly with the eventual formation of longer length-scale structures. This is 

160 shown in Fig. (6). The transport due to the electrostatic potential dominates substan- 

161 tially as depicted in the figure. Furthermore, in the steady-state, the nonlinearly coupled 

162 modes form stationary structures and D e ff saturates eventually. Thus, remarkably, an 
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163 enhanced cross-field transport level results primarily due to the emergence of large-scale 
W4 coherent structures in a turbulence dominated by the acoustic gravity waves. 

3. Concluding Remarks 

W5 In summary, we have investigated the turbulent properties of 2D nonlinearly interacting 

we AGWs by means of computer simulations. We found that the nonlinear equations gov- 

167 erning the dynamics of the fluid vorticity and density fluctuations in a nonuniform fluid 

us admit dual cascade, leading to a Kolmogorov-like energy spectrum. In such a cascade 

us process, we observe the formation of the short and large scale structures that are respon- 

™ sible for the cross-field (across the density gradient direction) diffusion of particles due to 

i7i random walk processes. The present results are essential for understanding the turbulent 

m properties of the atmospheres of the Earth and other planets. 
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Figure 1. Small scale fluctuations in the wave density fluctuations as a result of steady turbu- 
lence simulations of two-dimensional acoustic gravity waves. Forward cascades are responsible for 
the generation of small-scale fluctuations. The density perturbations are distributed isotropically 
in the two- dimensional domain. 

Figure 2. Evolution of potential fluctuations. Formation of large scale structures along the 
direction of the density inhomogeneity. 

Figure 3. Time evolution of the turbulent energy. The initial linear phase is followed by a 
nonlinear phase. The saturation follows at later times. 

Figure 4. Energy cascade rates. A rapid transfer of energy takes place during the early phase 
of the evolution followed by nonlinear saturation. In the latter, the energy cascade rates are 
saturated leading thereby to no net transfer amongst the turbulent eddies. 

Figure 5. The AGW turbulence exhibits a Kolmogorov-like k~ 5 ^ 3 power spectrum in the 
inertial range. 

Figure 6. Time evolution of the effective diffusion coefficient. 
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